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Non-equilibrium phenomena in the QCD phase transition 
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Physics Department, University of Athens, 15771 Athens, Greece 

Within the context of the linear a-model for two flavours, we investigate non-equilibrium phenom- 
ena that may occur during the QCD chiral phase transition in heavy-ion collisions. We assume that 
the chiral symmetry breaking is followed by a rapid quench so that the system falls out of thermal 
equilibrium. We study the mechanism for the amplification of the pion field during the oscillations 
of the (T-field towards and around its new minimum. We show that the pion spectrum develops a 
characteristic pronounced peak at low momenta. 



CT) • I. THE MODEL 

o : 

Experiments at RHIC and LHC are expected to probe many questions in strong interaction physics. One major 

area of interest concerns the chiral phase transition. For given baryon-number chemical potential fj, there exists a 

critical temperature T cr above which the system lies in the chirally symmetric state. As the temperature decreases 

i—i , below T cr the system moves into the chirally broken phase. It is believed that, for two flavours and zero quark masses, 

■ there is a lst-order phase transition line on the (T,/i) surface at large jj, pj- This line ends at a tri-critical point 

fNJ \ beyond which the phase transitions become 2nd order. The line of 2nd-order transitions ends on the fi = axis. In 

the case of non-zero quark masses, the lst-order line ends at a critical point, beyond which the 2nd-order transitions 

are replaced by analytical crossovers. 

This phase diagram has been discussed within various frameworks. Our interest lies in the study of possible 

non-equilibrium phenomena that may occur during the phase transition. In particular we would like to study the 

^ possibility that the system falls out of thermal equilibrium through rapid expansion. This is a realistic possibility in 

' the framework of heavy-ion collisions. 

The scenario we have in mind assumes an initial thermalization at a sufficiently high temperature for the system to 

C*~) ' move into the chirally symmetric phase. The subsequent fast expansion generates deviations from thermal equilibrium. 

We model this process by a quench during which the volume of the system increases instantaneously by a certain 

• factor, with the number densities of the various particles decreasing by the same factor. 

We consider only the two lightest flavours and neglect the effects of the strange quark. As an effective description 

JL . of the chiral theory we use the cr-model 0, Q . The Lagrangian density is 
Mh. 

<D . i 

£ = 2 ( V 5 "* + d^n) - V(a, t?) (1) 



> 

o 



with the potential 



- 1—1 

>< . 

V(a, t?) = ^(a 2 + 7T 2 - v 2 ) 2 + ^(cr 2 + n 2 - 2va + v 2 ) . (2) 

The last term in the potential accounts for the explicit chiral symmetry breaking by the quark masses. The scalar 
field a together with the pseudoscalar field if — (tt + , tt , 7r~) form a chiral field $ = (cr, ff). When the symmetry is 
restored at high temperatures in the absence of the second term in the potential, the system lies in the symmetric 
state (cr) = (ir) = 0. However, in the presence of the explicit symmetry breaking term in J5J, the expectation 
value of the c-field never vanishes, so that the chiral symmetry is never completely restored. At zero temperature 
and chemical potential the expectation values of the fields become (a) = v = f n and (7?) = 0, where f n is the 
pion decay constant. We fix the parameters of the Lagrangian using the phenomenological values ra,K 139 MeV, 
m a = y/2\ 2 f% + ml w 600 MeV, v w 87.4 MeV, which yield A 2 ~ 20. 
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II. EQUATIONS OF MOTION 

The equations of motion resulting from Q arc: 

<Vct + A 2 (ct 2 +k 2 -v 2 )a + ml<T = vm 2 

d^n + X 2 (a 2 +n 2 -V 2 )TT + mlTT = 0. 

We neglect the fluctuations of er, while we treat tt(x, t) as a quantum field: 



(3) 



(4) 



The creation and annihilation operators a\ v ,ak,v are defined in the interaction picture at the vacuum corresponding 

to the minimum v of the potential, and /&(t) are the mode functions of the pion field. 

We work within the framework of the Hartree approximation. This leads to the effective replacement Q 

• ir 2 (x,t)ir(x,t) -> §(t? 2 (X)) w(x,t), 

with (ir 2 (t)) — (ir 2 (t))/3 for i=l,2,3. The second approximation results from the replacement of the term nfirj by 
3(7r 2 )7Tj for j — i, and by (TT 2 )TTj for j ^ i. Substituting the above approximations into Q we get: 

a(t) + X 2 (a 2 (t) + {tt 2 (<)) - v 2 ) alt) + m 2 a(t) = vm\ (5) 



hit) 



k 2 -\ 2 v 2 + \ 2 a 2 (t) + ^\ 2 (jr 2 (t))+ml 



hit) = o. 



In 0,© (t? 2 ^)) is given by 



— ^/**(*)/*(*)- 



(2^ 



The particle density per momentum mode, for each component i of the pion field is : 



/ 



\fuW 



Ut) 



\ 



i 

2' 



with uij. = \Jk 2 + m* . For the total number of pions of all species (ir + , n° , 7r ) we have: 



E 



(6) 



(7) 



(8) 



V (2tt)3 

where is the volume of our system, that is the volume of the fireball in a heavy-ion collision experiment. 



(9) 



III. INITIAL CONDITIONS 

Our choice of the vacuum at v as our reference state has the advantage that the particle interpretation of the field 
w(x,t) is close to the experimentally observable quantities. It requires, however, some care with respect to our choice 
of initial conditions for the evolution described by equations (JHJ, ©■ We assume that the fireball created by the 
collision is initially in local thermodynamic equilibrium, or it has been separated in Disoriented Chiral Condensates 
(DCCs), each one in its own local thermodynamic equilibrium. If the second case is realised, our treatment applies 
to the interior of one of these DCCs. The initial expectation value <j\ that we consider for the c-field is small, but 
non-zero, because of the explicit chiral symmetry breaking. For our calculation we use <j\ = O.lu at T = 140 MeV. 
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For the pions we expect initially a thermalised gas that follows a Bose-Einstein distribution with 



eq 



1 



e t — 1 



(10) 



We assume the dispersion relation around the vacuum at a — v: 



even though the pion mass depends 



on the temperature. The justification for this approximation is provided by the explicit study in |7J of the effective 
pion mass during the process of chiral restoration. There, it is shown that stays approximately constant from 
T = up to T ~ 100 MeV. 

The mode functions fki(t), for a configuration corresponding to a non-interacting pion gas in thermal equilibrium, 
can be taken 



req _ 

Jki — 



,eq 



1/2 



— iu>kt 



(11) 



in agreement with (JHJ). In the following we restrict our analysis to large occupation numbers, and neglect the factor 
1/2 related to the zero point energy. 

In our simplified scenario we assume an instantaneous expansion of the fireball by a volume factor A (a quench). 
This means that the number densities of the pion gas must be reduced by the same factor. In addition, in order to 
be consistent with the conservation of energy, the initial value of the a-field has to change according to the relation 



V(a A ) 



A 



(12) 



where o\ is the value before the quench, and a\ the one after. We point out that this assumption is rather crude as 
it neglects possible fluctuations of a. However, it satisfies the minimal requirement of energy conservation. 
The above discussion implies that the physically motivated initial conditions for the evolution of the fields are 



fki(o) = \ -r 




/«(°) = V A 



(13) 



and 



(7(0) = <7A , <j(0) = 0. 



(14) 



These initial conditions differ from the ones assumed for particle production through inflaton decay in cosmology 
and in some wor ks on QCD HOI, 

as in those cases the initial particle number is taken to be zero. 



IV. NON-EQUILIBRIUM EVOLUTION 

Equations (JSJ) and JSJ, determining the evolution of the a and the pion field, constitute a non- linear integro- 
diffcrcntial system that is not solvable analytically. The main difficulty is related to the presence of the term (n 2 (t)) 
that mixes all the pion mode functions. If we ignore this term, it is possible to decouple the evolution of a from that 
of the pion field. We introduce the variable rj(t) = a(t)/v, and absorb X 2 v 2 in a new rescaled time variable r = Xvt. 
We also introduce the small parameter e = m 2 / (2X 2 v 2 ) <g; 1. In this way we obtain 

V"(r) + rf{T) + (2s - 1)t?(t) - 2e = 0, (15) 

where primes denote differentiation with respect to r. 

In zeroth order in e, i.e e = 0, this equation can be solved in terms of elliptic functions, with the result 

V(r) = y N , (16) 

dn^/T^M, 

where 



'l-7? 2 (0) 

1 v 2 (Q) • 



(17) 
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Here dn stands for the known Jacobi elliptic function and q is its modulus. (We do not use the standard symbol k in 
order to avoid confusion with momentum.) We point out that the solution (|16fl holds only for the range < 77(0) < 1, 
relevant to our case. For 77(0) > 1 the solution is expressed in terms of different elliptic functions. 
Inserting the solution Ijltjll into the equation for the pion mode functions, we obtain 



X 2 v 2 



ry 2 (0) 



dn 2 



1 - 



2 ' y 



2s- 1 



/w(r)=0. 



(18) 



This is the Lame equation, for particular values of the coefficients which make it solvable in terms of Jacobi functions 

MM- 

We first investigate the zeroth order solution (e = 0). As we show in the Appendix A, the solution /fci(r) of IjlSjl is 
a quasi-periodic function that can grow with time. Exponential amplification is obtained if the pion momentum lies 
in one of the zones: a) —00 < k 2 /(X 2 v 2 ) < 0, b) ry 2 (0) /2 < k 2 /(\ 2 v 2 ) < 1 — ?7 2 (0)/2. In our case the first zone is 
unphysical. Substituting a(t) = vrj(t), we obtain for the remaining amplification zone 



A 



A 



_ a ( )<fc<_V2 W 2-a 2 (0), 



(19) 



a result valid in zeroth order in e — m 2 r /(2X 2 v 2 ). 

Obtaining a full analytical solution valid in first order in e is very complicated. However, it is trivial to take into 
account the term 2e in equation l|18fl . Its consideration simply leads to the replacement of k 2 /{X 2 v 2 ) by 2e+k 2 / (X 2 v 2 ). 
The amplification zone becomes 



^a 2 (0)-m 2 <k< ] J^(2v 2 -a 2 (0)) 



(20) 



We emphasize that this result is not the complete answer in first order in e, as equation i|16fl is not an exact solution 
of (|15|) for £ ^ 0. However, the amplification zone given by equation l|2U|l agrees very well (within 4%) with the 
amplification observed through the numerical solution of the equations. 

If ct(0) is sufficiently small (er(0) < ^J^m^ ~ 0.5 v) the lower end of the amplification zone i|20|) is at k — 0. For 

larger <r(0) the zone becomes narrower. In the extreme case cr(0) ~ v (when the cx-field starts its evolution very close 
to its minimum), the zone shrinks to a point 



This is the pion momentum in the decay a — > 2tt. 

The form of equation ||SJ) for the pionic modes indicates the presence of two regimes in their evolution: 
a) Shortly after the quench the various modes have an effective mass term 



m lff 



k 2 - AV + AV(0) + ^A 2 (tt 2 (0)) + m 2 . 



(21) 



One expects the exponential growth of the low-momentum modes for which the mass term is negative. This phe- 
nomenon is characterized as spinodal decomposition 0, . 
b) At later stages the mass term becomes 



m lff 



(t) 



AV + AV(i) + ^A 2 (7f 2 (<))+m 2 



(22) 



with a a(t) a quasi-periodic function. If the term 



l (t)} is negligible the evolution equation becomes a Lame or 



Mathieu equation (see Appendix B for a detailed discussion of this point and a related controversy) . One expects a 
resonance band with exponential amplification of f^. This phenomenon is characterized as parametric resonance. It 
has been studied in detail in relation to the reheating of the universe through inflaton decay 0, IE IE • 

The analytical solution that we presented above, as well as the numerical analysis in the next section, demonstrate 
that there is no well-defined boundary separating periods in the time evolution during which one of the two mechanisms 
dominates. In general, there is a fast amplitude growth for the modes in a certain momentum zone (such as in equation 
d20l0 . However, the phenomenon can be a complicated convolution of both mechanisms we discussed above. 
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V. NUMERICAL RESULTS AND DISCUSSION 

We solve equations (J5J and |JB} numerically using a fourth-order Runge-Kutta algorithm for the differential equations 
and an 11-point Newton-Gotes integrator to compute the momentum integral: 

(7?2(<)> = / (0^ (<)/fe(<) = / ^ fc2 ^W/ fe W- ( 23 ) 

We calculate the pion density in 31? momentum space using (|SJ| and the total number of produced pions N to t using 
Furthermore, it is convenient to define the projected II? density p(k) through: 

/>oo 

N tot = / dkp(k). (24) 
Jo 

The pion gas is initially a fireball at thermal equilibrium with temperature T ~ 140 MeV and radius tq ~ 10 fm. 
The quench is described by an instantaneous expansion of this fireball by a factor A. We consider here three different 
values of the final radius corresponding to small (rj = llfm (A = 1.3)), intermediate (rf — 12.6 fm (A = 2)) and large 
(rf — 15fm (A = 3.375)) expansion of the initial fireball. The effect of the quench is incorporated in the normalization 
of the pion mode functions fk at t = 0. Performing the numerical integration we obtain the results presented in 
Figs. 1-8, which are close to our analytical calculations described in section IV. 

In Fig. 1 we depict the a evolution and the total pion number for expansion factor A = 3.375. The radius of the 
expanded fireball in this case is 15 fm. During the first two oscillations of a(t) the number of produced pions N to t 
increases fast, because of parametric resonance. This lasts about 6 fm and subsequently N tot just fluctuates around 
a mean value ~ 830. As our approach is semiclassical this solution describes the evolution of the system consistently 
only within the first stage (t ^ 6 fm). During this time the pion production can be described through the energy 
transfer from the oscillating a- field to the pionic field. In the following stages, our solution displays a periodic energy 
exchange between the two fields. A proper quantum treatment is necessary in order to describe the complete decay of 
the cr-field into free pions. However, the bulk of pions is produced during the first stage. For this reason we estimate 
the various observables (such as pion number and distribution) at a time t ~ 6 fm. 

In Fig. 2 the evolution of the distribution of pions p{k) for various times is shown. We observe a large enhancement 
of the spectrum at low momenta, and the formation of a zone with a peak at a specific fc-value (k max ) and a width 
characteristic of the non-equilibrium amplification. Comparing with Fig. 1 we observe that the energy transfer to 
the pions is translated to a decrease of the cr-oscillation amplitude. In Fig. 2 we also observe a slight sinking of the 
spectrum at large momenta, which implies an additional energy transfer from hard to soft pion modes through the 
mode- mode coupling (the (7? 2 (i)) term) in (6), effect that cannot be described within our approximations. 

This behaviour can be described very well using the analytical results of section IV. For an expansion factor 
A = 3.375 (and more generally for A > 3) the initial values fki(0) (according to (fH5|) 1 are very small so that 
(if 2 (0)) <C v 2 , justifying the conditions for parametric resonance: while a oscillates around its new minimum after the 
quench, it leads the pion mode functions in the momentum zone <|20ll to exponential increase. The numerical values for 
the limits of the momentum zone, whithin which the exponential amplification occurs, agree with a 4% accuracy with 
the analytical estimate of equation (|20|l . This is explained by the fact that {ir 2 (t)) is negligible initially in this case, 
so that the analytical solution gives a good approximation. As time passes, (ir 2 (t)} increases, so that eventually the 
parametric resonance ends. Even though this effect cannot be described within the analytical approach, its influence 
on the location of the amplification zone is very small. 

We also point out that the initial conditions for a(t) are determined by equation (|12|) . so that a large A corresponds 
to cr(0) closer to v. This large initial value of a shifts the pion effective mass at t — to larger values, leading to a 
strong suppression of spinodal decomposition in favour of parametric resonance. 

In Fig. 3 we plot u(t) and N tot (t) for A = 1.3. We observe that the pion production completes during the rolling 
down of the cr-field towards the new minimum. This process lasts for ~ 2 fm, and subsequently N to t oscillates in a 
manner similar to that in the case A = 3.375. Observables are now considered for times t ~ 2 fm. For small expansion 
factors (A < 1.5) the dominant amplification mechanism is the spinodal decomposition: during the first rolling down 
of the cr-field towards its minimum, the effective squared pion mass <|21ll is negative and the pion modes with small 
momenta are amplified. Furthermore, because of the large value of (tt 2 (0)}, the requirements of parametric resonance 
are not fulfilled in the second stage when cr(i) oscillates around the new minimum. As a result the amplification 
through parametric resonance is suppressed. We point out that the analytical solution of the previous section does 
not give a good quantitative approximation in this case, as the effect of the term ~ (Tr 2 (t)) cannot be neglected. 

For completeness, we show in Fig. 4 the evolution of the pion distribution for A = 1.3. Similarly to the previous 
case (Fig. 2) we observe an enhancement of the pion spectrum at low momenta. The lower limit of the ampilified 
zone is at k = 0, while the upper limit is time dependent. The peak of the particle density per momentum mode is 
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FIG. 1: The evolution of the <r-field as well as the total pion number Ntot(t) (scaled by a factor of 3) for initial radius 10 fm 
and final radius 15 fm (A = 3.375). 




FIG. 2: The density p(k) defined in equation il'J 1 fi using A = 3.375 at three successive times. 
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FIG. 3: The evolution of the cr-field as well as the total pion number Ntot(t) (scaled by a factor of 3) for expansion factor 



at k = 0. However this feature is modified at the level of one-particle momentum density <|24[) . because of the phase 
space factor Airk 2 that shifts the maximum to a non-zero momentum value. 

Lastly, for intermediate expansion factors (1.5 < A < 3) both amplification mechanisms are operative. During the 
first CT-rollmg, m e ff(t) is negative for low momenta and spinodal decomposition takes place. Subsequently, while a 
oscillates around its new minimum, (7? 2 (i)) is not very large, allowing for parametric resonance to occur as well. This 
behaviour can be clearly seen in Fig. 5: the pion production happens during both phases. At late times (7r 2 (i)) grows 
significantly, parametric resonance ends, and pion production completes. The time when this happens is ~ 4 fm. 

In this intermediate case both limits of the amplification zone are time dependent. As a result, the pion momentum 
spectrum changes with time. In fact, we get an oscillating zone in momentum space within which the amplification 
takes place. A simple explanation of this behaviour can be obtained if we observe that the addition of a constant 
term within the brackets in equation l |18(l shifts the amplification zone given by equation (|20(1 , to the left if the added 
term is positive, or to the right if it is negative. The addition of a fluctuating term ((n 2 (t)}) is expected to force the 
edges of the amplification zone to oscillate. It must be pointed out that, after the inclusion of the back-reaction term 
(ir 2 (t)) in the equation for a, the solution (|16|) is no longer valid. However, in this intermediate case the modification 
of the solution is not expected to be dramatic. The expected behaviour is observed in Fig. 6. The pion spectrum is 
enhanced and its maximum is shifted towards lower momenta as in the previous two cases. However, the location of 
the maximum of the spectrum and the limits of the amplification zone oscillate. 

A common property of the momentum spectrum of the produced pions, for all the expansion factors we considered, 
is the appearance of an amplification zone in the region of low momenta. This results in the production of a relatively 
large number of soft pions during the evolution of the a-w system. As we have already mentioned, our semiclassical 
approach is not reliable for long times. We estimate the number of produced pions through the pion distribution at 
the beginning of the stage in which our solution predicts a periodic interchange of energy between the cr-field and 
pion sectors: ~ 6 fm for A = 3.375, ~ 4 fm for A = 2 and ~ 2 fm for A = 1.3. The effect of the amplification is quite 
strong, and persists even if we perform an averaging of the solution over times larger than the ones listed above. The 
energy density p(k) is clearly distinguishable from that of a thermal system because of the presence of a narrow peak 
at low momenta. 

The formation of the momentum amplification zone reflects a clustering of the pion momenta in the low fc-region. 
The dynamical character of this cluster, which leads to the presence of a peak at a definite momentum value, ac- 
companied by a width of typical size, in the projected one-particle density, is crucial. This characteristic form can 



A = 1.3. 
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FIG. 4: The density p(k) using A = 1.3 at three successive times. 




— i 1 1 1 1 1 1 1 1 1 

2 4 6 8 

t(fm) 

FIG. 5: The function a(t) as well as the total pion number N tot (t) (scaled by a factor of 3) for an expansion factor A = 2. 
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FIG. 6: The density p(k) for the case of intermediate expansion (A = 2) at four successive times. 



be directly observed in an inclusive analysis of the pion momenta produced for example in A + A ultrarelativistic 
collisions. The geometry of this cluster is expected to be rather simple: no substructures at smaller scales are present. 
It is worth comparing the situation here with the corresponding formation of particle clusters in a critical system. 
In the latter case the clusters have statistical character. Therefore, at the level of one-particle inclusive density in 
momentum space no peaks (beyond the kinematic ones) appear. The geometry of these critical clusters is that of a 
random fractal: self-similar structures at different scales occur, leading to a characteristic power-law dependence of 
the factorial moments on the resolution scale M [l7| . 

In order to explore the resonant pion production at a phenomenological level, we have generated, through a Monte- 
Carlo simulation of equations (fO^) . H24JI . data sets consisting of a large number of events. For each value of the expansion 
factor A that we considered, we produce 1000 events. We then calculated the corresponding second factorial moment 
in transverse momentum space as a function of the resolution scale M [l^. O ur results are given in Fig. 7. We observe 
a conventional behaviour of the moments (saturation at high resolution scales) reflecting the absence of self-similar 
structure at different scales in the formed pionic clusters, for all expansion factors. 

One could also consider the inclusive density of the transverse momentum for each set of 1000 events. We have 
calculated this distribution for the case A = 3.375. The result is displayed in Fig. 8, where we clearly observe the effect 
of the amplification zone. In the same plot we show the experimentally observed transverse momentum distribution 
of the charged pions produced in Pb + Pb collisions at 158 GeV/n (NA49-SPS) |l9|. The experimental results are 
very close to a thermal system in transverse momentum space. The large deviation of the transverse momentum 
distribution of the pions produced through the amplification mechanisms described in the present work, from the 
corresponding distribution of a conventional thermal system, is an experimentally accessible signal for the appearance 
of non-equlibrium phenomena, associated with chiral phase transition, in the pion production during a heavy-ion 
collision. 



VI. CONCLUSIONS 



In this work we studied non-equlibrium phenomena related to the chiral QCD phase transition. We found that, 
within the scenario of quench, the dynamics of the system may amplify pion modes in a certain momentum zone 
while the cr-field moves towards and around the (chirally broken) vacuum. These phenomena produce a significant 
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FIG. 7: The second factorial moment for the three different cases of large (A = 3.375), intermediate (A = 2) and small 
(A = 1.3) expansion. For each case we have analysed 1000 Monte-Carlo events generated according to equations 1241 1. 
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FIG. 8: The transverse momentum distribution of the finally produced pions in the case of large expansion (A = 3.375). For 
comparison the fcT-distribution for a system at thermal equilibrium (Pb + Pb at 158 GeV/n) measured at the NA49 experiment 
(CERN-SPS) is also displayed. Both distributions are normalized to the same total pion multiplicity. 
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enhancement of the pion spectrum at low momenta. In our approach we took account of two main features of the 
physical system: the explicit symmetry breaking through the non-zero pion mass, and the local thermal equilibrium 
for the initial pion gas. In this sense our study is close to the real world. As a consequence, the pion production in our 
model is not as pronounced as in other models [^.l^llOj. However, the momentum distribution of the produced pions 
possesses fingerprints of the zone structure of the amplification regions in momentum space: a cluster of pions with 
low momenta is formed. Although visible at the level of single particle density, this cluster does not have substructures 
at different scales. It is, therefore, distinguishable from particle clusters created by different mechanisms (equilibrium 
phase transition). Furthermore, the transverse momentum spectrum acquires a specific non-thermal peak that could 
be a clear signature of these non-equilibrium phenomena in a collision experiment with relativistic heavy ions. 

It would be interesting to examine the modifications in the scenario we studied, if the finite time neccessary for 
the quench is taken into account. Additionally, one should also consider the consequences of the cylindrical form of 
the expanding fireball for the observed pion spectrum. Such a geometry is closer to the conditions occuring in A + A 
experiments at very high energies. Another important aspect is the case when the initial state is critical, i.e. large 
self-similar fluctuations of the er-field are present. In this case, the interesting problem is to explore the evolution of 
the critical state towards the freeze-out phase. The study of these additional issues is left for future investigations. 

Appendix A 

We present here the derivation of some of the results used in the main text. Equation l|18fl with e = has the form 

k 2 



with 



X 2 v 



+ 7?»-l 



/fci(r) = (25) 



1 - 



V 2 (0) 



(27) 



Here t](t) is a doubly periodic function with periods 2uj = 2K(q) / \J 1 — (real) and 4u>' = AiK'(q)/ \J 1 — 
(complex), where K(q) is the complete elliptic integral of first kind and K'(q) is its complementary. 

Equation l|25|) is the Lame equation, for a particular value of the coefficients that make it solvable in terms of Jacobi 
functions. It is convenient to express the Jacobi function dn in terms of the Weierstrass V function through |12| 



dn (u,q) 



1 (u+K (q)+iK T (qY) 

l~q 2 



ex - e 2 



u + K(q) + iK'(q) 



(28) 



■ The parameters ex, e%, e$ are the roots of a cubic equation whose coefficients depend on the two 



where q = 

periods of the Weierstrass V function 0. They obey the relation ex 
choose ex — e$ = 1 — v ^ . Thus we obtain 

3ei 



i, 

e-i + e 3 = 0. Without loss of generality we 



3e2 



1 



3„2 



n 2 (o) 



3e 3 = §77 2 (0)- 2. 



Inserting (|28|l into (|25|l and using the definitions of the periods lu,lu' 

fki(r) + 

The solutions of (|30|l are of the form ^lj 



we find 
fki(r) = 0. 



, = g^T+^W + w(fc)) e _ TC(w(k)); 

a(j + uj + uj') 



(29) 



(30) 



(31) 
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where <j(x) and £(x) are the known Weierstrass functions defined as 



u 

and 



cr(u) = uexp | y 



C(z)--)rfzf., (88) 



and the complex number w(/c) is defined through the transcendental equation 

*W~3& + S- (34> 

Because of the quasi-periodicity property of the a-function 

tr(a: + 2w) = e 2 (^+")C("), 

we obtain from l(3"T|) 

f ki (r + 2(j) = /fci ( T ) e 2[w(fc)CH— C(w(fc))]_ (35) 
As a result, for w(fc) for which the exponent 

At(w(fc)) = 2 [w(fc)C(w) - wC(w(A))] (36) 

becomes real, we have an exponential amplification of /fci(r). For values of w(fc) which lead to imaginary /x(w(fc)) we 
obtain oscillatory behaviour. 

To determine the corresponding /i(w(fc)) intervals one has to consider the equation l|34|l that maps the real k 2 / (X 2 v 2 ) 
axis onto the sides of the fundamental square in the w plane [Ed 

• If -oo < ^ < 0, then w(fc) = £ R with < (3 < oj. (A) 

• If < < 2^21 , then w(fc) = uj + ia (a e R) with < a < -id. (B) 
. if A < < i _ 2fM f then w[k) = d + i/3 with w > (3 > 0. (C) 

• If 1 - 2!M < ^ < oo, then w(fc) = with -id > a > 0. (D) 

The appearance of 77(0) in these four intervals results from the dependence of the roots ei, e2, e3 in l)29fl on ?7(0). Using 
the properties of the £- function [12| we can easily derive that /x(w(fc)) has a non-zero real positive part for the cases 
(A) and (C), while it is purely imaginary for (B) and (D). 

Appendix B 

We derive the Mathieu equation starting from the Lame equation. We use the expansion d] 

1 7T 



dn(u, g ) 2q'K(q) 

where 



(37) 



Q = e *■<«) 



is the elliptic nomc. This expansion holds for Im ( 2K\q) ) ^ s 71 "! 111 (t^t) ^3" Q can ^ e ex P an ded as 



Q^e-^^W- = A + 2A 5 + 15A 9 + 150A 13 + • • • (38) 
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where 



A = 



K(q) can be also expanded as 



1 l-y/g 7 = 1 l-(l-g 2 ) 1/4 

2 1 + y/? ~ 2 l + (l_ g 2)l/4- 



2™n! 



(39) 



Keeping only the first terms in the expansions we obtain 

V 2 (0) V 2 (0)n 



ff(r) = 



2(1 -w-y*) 



1 - 8(A + 2A 5 ) cos 



2tJ1 



n 2 (0) 



i + k 2 + ^ 4 



(40) 



with q 



1 ^(oj 1 ■ Substituting this in the Lame equation 1251) we find 



r ^2 



dz 2 



with 



A 
D 



A - 2Scos2z 



2s- 1 



/fci(2) = 



77 2 (0)tt 



(41) 



= 2 



X 2 v 2 ' ~ * ' 2(1 

77 2 (0)7T 



h 2 - ^) 



(i 



3^ 



(A + 2A S ) 



n 2 (0) 
2 



1«2 



64 



(42) 



Equation (|41l) is the Mathieu equation. It is a good approximation of the Lame equation for Q < 1, i.e. for q -C 1, 
since the elliptic functions approach the usual trigonometric ones for q — > 0. However, in the case Q ~ 1, i.e for 
q ~ 1, all higher-order terms in (|37J) must be taken into account. The simple approximation that leads to the 
Mathieu equation is not valid. There is some controversy in the literature [l6l | on whether the Mathieu equation gives 
a reliable approximation in the studies of parametric resonance during inflaton decay. 
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